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The roton in a Bose-Einstein condensate is computed, both near and far from a Feshbach res- 
onance. A low-density approximation is made, allowing for an analytic result. A Monte Carlo 
calculation shows that the roton is larger than predicted by the low-density approximation, for the 
upper range of densities considered here. The low-density approximation is applied to superfluid 
4 He, roughly reproducing the results of previous Monte Carlo calculations. 

PACS numbers: 



A roton, first predicted by Landau [l| for superfluid 
4 He, is an excitation in a Bose-Einstein condensed fluid, 
characterized by a minimum in the excitation spectrum 
uj(k). By Feynman's relation 0, this minimum also cor- 
responds to a maximum in the static structure factor 
S(k). This peak in S(k) occurs for excitations whose 
wavelength 2w/k is equal to the characteristic wavelength 
of density fluctuations in the ground-state wave function 
of the quantum fluid. The peak in S(k) exceeds unity, 
the value for an uncorrelated gas. 

Feynman described several possible microscopic pic- 
tures of a roton in superfluid 4 He. He found that the 
most likely description is that the roton is analogous to 
a single atom moving through the condensate, with wave 
number k close to 27r/n -1 / 3 , where n -1 / 3 is the mean 
atomic spacing. 

The roton in superfluid 4 He was calculated by a Monte 
Carlo technique 0, 0] using a many-body Jastrow wave 
function [j| , which is of the form 

N 

i>= n f(\Ti- rj \), (i) 

j>i=l 

where the pair function f(r) should be determined for the 
quantum fluid under consideration. The wave function 
l|T| deviates significantly from unity if any two atoms 
become very close to one another. For superfluid 4 He, 
the two-particle correlation function g(r) for this wave 
function has fluctuations with a preferential length scale 
of n -1 / 3 , which results in a peak greater than unity in 
S(k) near k f» 27r/n~ 1 / 3 , in agreement with Feynman's 
picture. 

Both Feynman's result and the results of the Jastrow 
wave function roughly agree with measurements [6( of 
S(k) and u{k). 

The Jastrow wave function 0) has also been consid- 
ered for a low-density Bose-Einstein condensate (BEC) 
, and has been used to compute various properties of a 
high-density BEC 0. We employ a Jastrow wave func- 
tion to compute the roton for a low-density BEC. By an 
analytic calculation, we find that the hard-sphere size r Q 
of the atoms determines the location of the roton, rather 
than n -1 / 3 . We see that Feynman's view that n -1 / 3 



is the relevant length scale is the special case of high- 
density, for which r a w n -1 / 3 . 

In a low-density BEC, Feynman's view of the roton as 
a single atom moving through the condensate seems nat- 
ural. In general, the excitation spectrum of a low-density 
BEC is of the Bogoliubov form 0, [j| , which consists of 
phonons and single-particle excitations. We find that the 
roton occurs on the single-particle part of the spectrum, 
at A; ~ 8/7ra, where a is the s-wave scattering length. 

Although the general form of the wave function (JJJ can 
be applied to both superfluid 4 He and BEC, the wave 
function is fundamentally different for these two quan- 
tum fluids. For both of these fluids, below some temper- 
ature T s - wave , the thermal de-Broglie wavelength X^b 
becomes longer than the characteristic length scale R of 
the interparticle potential. Below T s — wave , all scattering 
processes except for s-wave scattering become negligible. 
Since a BEC is a dilute gas, n -1 / 3 , which is typically 
about 1500A, is much larger than R. Therefore, the crit- 
ical temperature T c for quantum degeneracy, at which 
XdB becomes comparable to n -1 / 3 , is much lower than 
Ts-wave- Therefore, only s-wave scattering plays a role 
here for a BEC. In contrast to BEC, superfluid 4 He is 
a relatively dense liquid, for which n -1 / 3 is comparable 
to R. Therefore, the temperature T s - wave also marks 
the transition T c to quantum degeneracy. Even below T c 
for superfluid 4 He, partial waves other than the s-wave 
contribute to the wave function. 

For many properties of a BEC, a is sufficient to de- 
scribe the interparticle potential, and the form of the 
potential does not play a role [lfl Hi] ]. To quantita- 
tively describe the roton however, the s-wave scattering 
wave function must be known to interparticle separations 
somewhat smaller than R, where the details of the po- 
tential are relevant. 

Rotons have been pred icted in a BEC in the presence 
of optical fields |~L2l Il3|. as well as in a dipolar BEC 
(l4| . In this work, we find that a roton occurs in an 
unperturbed BEC without any dipole interaction. We 
consider the enhancement of the roton near a Feshbach 
resonance 0, [l^, but the physics is qualitatively the 
same as in the unperturbed case. 

When we refer to a roton, we are referring to a peak in 
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S(k), rather than a minimum in ui(k). By the Feynman 
relation, the peak in S(k) computed here for a BEC is 
not steep enough to produce a minimum in u>(k), for the 
range of densities considered here. 

We use a low-density approximation |5j to compute 
g(r) and the roton for a BEC with a positive scattering 
length. The low-density approximation is compared to a 
Monte Carlo calculation. 

For a BEC of alkali atoms, the potential can be taken 
as the van der Waals potential —Cs/r 6 for atomic sepa- 
rations r greater than a few Angstroms 0, 0] . For this 
potential, R is given by (mC^/H 2 ) 1 / 4 , which ranges from 
about 50 to 100A . The wave function for s-wave 

scattering between two particles in the limit of zero en- 
ergy is given by ipQ (r) , which is the solution of the radial 
Schrodinger equation with the van der Walls potential. 
ipe(r) can be written in terms of elliptic integrals [l7| . 
For r > R, this wave function is very close in form to 
1 - a/r [lfl,|2I|. For the case of a > R. ibp(r) therefore 
obtains large negative values for r -C a |2lj . 

We cannot use tp&{r) as f(r) in Q, because even for 
r 3> n -1 / 3 , ipQ (r) is significantly less than unity. This is 
non-physical, in the sense that no matter how large the 
gas for fixed density, the value of the many-body wave 
function (JJJ depends on the size of the gas. To account 
for many-body effects, we use the following pair function 
which goes to unity for r > n -1 / 3 0>E|- 

/(r) = (iWMM 3 ) (r<n-V3) 

To compute (0), R is needed. Throughout this work, we 
use R = 0.05 n- 1 / 3 . This is a typical experimental value, 
and the results here are rather insensitive to R. 

Equation JJJ with J5J is shown schematically in Fig. 1. 
To aid in visualization, Fig. 1 shows the wave function 
squared in one dimension, as a function of the position 
x\ of atom number 1. The positions of all of the other 
atoms, such as atoms b through e, are fixed. As long 
as atom number 1 is far from the other atoms, the wave 
function has the constant value %j) Q . This value is deter- 
mined by the positions of the atoms other than X\. The 
rapid oscillations in ip§(r) for r < R appear in Fig. 1 as 
dark vertical bands, when atom number 1 is very close 
to another atom. 

The correlation function g(r) gives the unconditional 
probability of two atoms being at a distance r. g(r) is 
related to the pair function f(r) by [j| 

N 

Jdr 3 ...dr N n / a (k<-'il) 

g(\r 1 -r 2 \) = V 2 ^ , (3) 

fdr^drv J] f 2 (\n-rj\) 

j>i=l 

where the integrals are over the volume V. © can be 
evaluated starting with the first integral in the denomi- 
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nator, J dr± FT / 2 (|i"i — r ^ | ) . This can be visualized as 

the integral over the function shown in Fig. 1 for one di- 
mension. Neglecting three-body interactions greatly sim- 
plifies this integral. Three-body interactions are rare for 
small values of the gas parameter na 3 , assuming that 
the range of three-body interactions is of the same order 
of magnitude as the range of two-body interactions. A 
three-body interaction is represented in Fig. 1 by the 
points x c and Xd, where atoms 1, c, and d interact. 
Neglecting such interactions, the integral is a function 
of the volume v indicated by the shaded region in Fig. 
1. The integral is then given by V — (N — l)v, where 
v = V — J drf 2 (r). This result is independent of the 
positions of the . 




FIG. 1: Schematic one-dimensional representation of the Jas- 
trow wave function. The dependence on the single dimension 
Xi is shown. The overall scaling of the curve is given by the 
factors in the wave function not involving xi. The labeled 
values of x\ correspond to the positions of atoms other than 
atom 1. x\ varies over the volume V. 

Evaluating all of the integrals in © similarly to the 
first yields g(\ ri - r 2 |) » V 2 / 2 (|ri - r 2 \)[(V - v)V]-\ 
Since V ^> v, we obtain the result of the low-density 
approximation || 

g(r) « f(r). (4) 

The result Q with J3J is indicated by the solid curve 
in Fig. 2 for na 3 = 2.2 x 10~ 4 . This value of na 3 is 
an order of magnitude greater than typical experimen- 
tal values without a Feshbach resonance. The result for 
na 3 = 0.01 1 is also shown in the figure. 

While @ is useful for obtaining analytic results, a 
more accurate computation can be made by the Monte 
Carlo technique described in Ref. |3| . This technique ef- 
fectively evaluates @ by using a Metropolis algorithm to 
randomly choose configurations of the N atoms (N =100 
here), according to the probability distribution given by 
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||TJ with (J2J , and computing the distribution of distances 
between the atoms, with periodic boundary conditions. 
This distribution, averaged over many likely configura- 
tions, is proportional to r 2 g(r). The result of the Monte 
Carlo computation is shown in Fig. 2 for na 3 = 2.2 x 10 -4 
and 0.011. These results were obtained with 9 x 10 8 and 
2 x 10 7 iterations, respectively. 
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FIG. 2: The two-particle correlation function for a BEC. The 
solid and dashed curves are the low-density approximation Q 
with with na 3 = 2.2 x 10~ 4 and 0.011, respectively. The 
dotted curve is the Monte Carlo result for na 3 = 0.011. For 
na 3 — 2.2 x 10~ 4 , the Monte Carlo result is indistinguishable 
from the solid curve. 
0.011 only. 



The inset shows the curves for na 



The small oscillations for small r shown in Fig. 2 are 
negligible in the computation of the static structure fac- 
tor. These are the oscillations for r values less than that 
of the first large peak below r/a — 1, in the solid and 
dashed curves. To save computing time these oscillations 
are not included in f 2 (r) in the Monte Carlo computa- 
tion. 

For a > 1.3R (corresponding to na 3 > 2.7x 10~ 4 in this 
work), The peak in / 2 (r) located at r < R is greater than 
unity. The dashed line of the inset of Fig. 2 shows an 
example of such a peak, for a = 4AR. Because this peak 
is greater than unity, a cluster of atoms is the most likely 
configuration. Since we are interested in the un-clustered 
phase, we do not let the Monte Carlo computation pro- 
ceed long enough for the clear transition to clusters to 
occur. 

The static structure factor S(k) is given by 1 + 
n J[g(r) — l]e lk ' r dr, which in general can be written 



S(k) = l + 4irn / drr 2 [g{r) - 1] 



sin(kr) 



ki 



(5) 



We use (JSJ) to compute S(k) for the low-density approx- 
imation (0J with J2J), as indicated in Fig. 3. The height 



S(k r ) and location k r of the roton are indicated by the 
solid curves of Fig. 4. S{k) in Fig. 3 has several maxima, 
the tallest of which is taken as the roton. As na 3 is var- 
ied, the peak which is taken as the roton varies, resulting 
in the jagged appearance of the solid curves of Fig. 4. 

For the Monte Carlo computation, S(k) is found by 
inserting g(r) such as is shown in Fig. 2 into JSJ. The 
results are indicated in Fig. 3 by the dash dotted and 
solid curves. For small na 3 , the height and location of the 
roton are seen in Fig. 4 to be the same for the low-density 
approximation and the Monte Carlo result. The low- 
density approximation is therefore valid for small na 3 . 
For the larger values of na 3 , the height of the roton in 
the Monte Carlo calculation is larger than that of the 
low-density approximation, as seen in Fig. 4a. 

By applying the low-density approximation to the limit 
of a 3> R, we can obtain an analytic expression for S(k) 
for R -C a <C n^ 1 / 3 . In this range, f(r) can be taken as 
1 - a/r. By (@J and ©, 



S(k) = 1 + 4irna 3 [Tr(kay 1 /2 - 2(ka)~ 



(6) 



This small-R limit is indicated by the dotted curve in Fig. 
3. The height of the roton in this limit is S(k r ) = 1 + 



■K 3 na 3 /8, and the location is given by k r 



These 



values are indicated by the dashed curves of Fig. 4, which 
as expected, agree with the low-density approximation 
(solid curves) for a ^> R, where nR 3 is indicated by the 
dotted lines. As seen in the above expression for k r , the 
location of the roton is determined by o, rather than by 

n -l/3 
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FIG. 3: The roton in a BEC, computed by 0. The dash 
dotted and solid curves are the Monte Carlo results for na 3 = 
0.011 and 2.2 x 10 -4 respectively. The result of the low-density 
approximation Q with © for na 3 = 0.011 is indicated by the 
dashed curve. The result of the low-density approximation 
for na 3 = 2.2 x 10 -4 is indistinguishable from the solid curve. 
The dotted curve is the small-R limit, given by JSJ. 



For na 3 — 0.011, the height of the roton given by the 
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Monte Carlo calculation is S(k r ) — 1.08, as shown in 
Figs. 3 and 4a. Measuring this 8% effect could be exper- 
imentally feasible, na 3 — 0.011 could be attained by a 
Feshbach resonance [2(|. It should be noted though, that 
for this relatively large value of na 3 , the form of (J2J is 
only approximate, so the expressions for the height and 
location of the roton should be considered as estimates 
only. 

Due to phonons, the true S(k) is proportional to k for 
k < where f — 1 = a~ 1 \f&ima? is the inverse healing 
length. The curves shown in Fig. 3 do not show this 
linear behavior for small k because the wave function JJJ 
with Q does not have the long-range correlations of a 
phonon 0,0]. 

For superfluid 4 He the three inverse length scales 
27r/n -1 / 3 , k r , and are roughly equal. This is also 
true for a strong roton in a BEC. As na 3 increases, both 
k r and £ _1 approach 2n/n~ 1 / 3 . More precisely, the ra- 
tios of and k r to 27r/n~ 1 / 3 are ,/2/ir(iia 3 ) 1 / 6 and 
approximately 4/7r 2 (na 3 )~ 1 / 3 , respectively. The latter 
ratio implies that an appropriate measurement system for 
measuring a roton should be able to probe wavelengths 
somewhat shorter than n -1 / 3 . 
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FIG. 4: The height (a) and location (b) of the roton in a 
BEC, as a function of na . The circles are the Monte Carlo 
result. The solid curves are the low-density approximation. 
The dashed curves are the analytic small-R limit. The dotted 
lines indicate ni? 3 . 

As a demonstration of the validity of the low-density 
approximation Q), we use (J3J to compute various prop- 
erties of superfluid 4 He found in 0] and 0] by means of 
Monte Carlo calculations. For superfluid 4 Hc, nb 3 sa 0.4 
where b is the hard-sphere size. Therefore, this is a strin- 



gent test of the low-density approximation. We find that 
the location of the roton for superfluid He is given by 
k r = 5.3/6, which differs from the Monte Carlo result 
by 6% 0, EJ. We find that the height of the roton for 
superfluid 4 He is S(k r ) fall 0.3n& 3 , which gives a roton 
of S(k r ) = 1.1, compared to S(k r ) — 1.2 from the Monte 
Carlo calculation HQ, and S(k r ) = 1.5 from experiment 
6] . The order-of-magnitude agreement between the low- 
density approximation and the Monte Carlo result for 
superfluid 4 He suggests that the low-density approxima- 
tion preserves the essence of J3Jl. 

In conclusion, we find the height S(k r ) and location 
k r of a roton in a BEC, for a range of densities. A low- 
density approximation is compared to a Monte Carlo cal- 
culation. The values of S(k r ) and k r given by the two 
methods agree for the lowest densities. For higher den- 
sities, the Monte Carlo calculation predicts an enhance- 
ment in S(k r ) of almost a factor of 2 over the low-density 
approximation. 

In contrast to the Monte Carlo calculation for super- 
fluid 4 He, the small-R limit gives explicit expressions for 
the height and location of the roton. 
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